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ABSTRACT 

We used the OGLE-II and OGLE-III photometry of red giants in the Large Magellanic 
Cloud to select and study objects revealing ellipsoidal variability. We detected 1546 candidates 
for long period ellipsoidal variables and 121 eclipsing binary systems with clear ellipsoidal mod- 
ulation. The ellipsoidal red giants follow a period-luminosity (PL) relationship (sequence E), 
and the scatter of the relation is correlated with the amplitude of variability: the larger the 
amplitude, the smaller the scatter. 

We note that some of the ellipsoidal candidates exhibit simultaneously OGLE Small Am- 
plitude Red Giants pulsations. Thus, in some cases the Long Secondary Period (LSP) phe- 
nomenon can be explained by the ellipsoidal modulation. 

We also select about 1600 red giants with distinct LSP, which are not ellipsoidal variables. 
We discover that besides the sequence D in the PL diagram known before, the LSP giants form 
additional less numerous sequence for longer periods. We notice that the PL sequence of the 
ellipsoidal candidates is a direct continuation of the LSP sequence toward fainter stars, what 
might suggest that the LSP phenomenon is related to binarity but there are strong arguments 
against such a possibility. 

About 10% of the presented light curves reveal clear deformation by the eccentricity of 
the system orbits. The largest estimated eccentricity in our sample is about 0.4. 

All presented data, including individual B observations and finding charts are available 
from the OGLE Internet archive. 

binaries: close - binaries: eclipsing - Stars: late-type - Magellanic Clouds 



* Based on observations obtained with the 1.3 m Warsaw telescope at the Las Campanas 
Observatory of the Carnegie Institution of Washington. 
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1. Introduction 

Ellipsoidal variability is a phenomenon often observed in close binary systems. It 
is caused by aspect changes of the component deviated from spherical symmetry 
by tidal interactions with its companion. A typical ellipsoidal light curve is 
sinusoidal-shaped, with two maxima and two minima per orbital period. Light 
variations should peak at orbital phases 0.25 and 0.75, when maximal area of 
the Roche-lobe is projected. In general the depths of the minima are not equal, 
which is the result of next order effects of the tidal distortion. 

Ellipsoidal variables are relatively rarely studied objects, and long period 
ellipsoidal variables are almost unknown. Morris (1985) gathered all ellipsoidal 
variables known at that time, and listed 20 confirmed and 20 suspected ob- 
jects of this type. Only one variable star from this catalog (T CrB) had period 
longer than 100 days. Since then, the number of known ellipsoidal variables 
has increased, but still, the list of confirmed ellipsoidal red giants is very short. 
However, the analysis of the ellipsoidal variations is a useful tool for study- 
ing various types of close binary systems: X-ray binaries, cataclysmic binaries, 
symbiotic stars, early type near contact systems. 

In recent years gravitational microlensing surveys have revolutionized our 
knowledge about variable red giants. Wood et al. (1999) presented period- 
luminosity (PL) diagram of the long period variables which showed series of 
distinct parallel sequences, marked with letters A-E. 

In the previous paper (Soszyhski et al. 2004), we selected and studied about 
15 000 red giants occupying most often the sequence A. We showed that these 
objects, named by Wray, Eyer and Paczynski (2004) OGLE Small Amplitude 
Red Giants (OSARG), constitute a separate type of pulsating giants, different 
than "classical" semi-regular variables (SRV) . 

In this paper we present a sample of variable red giants constituting PL 
sequence E, spreading below the sequences A-D. Wood et al. (1999) suggested 
that most of these stars are contact binaries. We studied photometry of these 
objects and noticed that, indeed, the majority of the light curves are probably 
ellipsoidal or eclipsing close binary systems. Our sample is, then, a natural 
extension of the catalog of eclipsing binary systems in the LMC prepared by 
Wyrzykowski et al. (2003). 

2. Observations and Data Reductions 

Observations presented in this paper were carried out with the 1.3-m Warsaw 
telescope at the Las Campanas Observatory, Chile, operated by the Carnegie 
Institution of Washington. Central parts of the LMC were added to the list of 
regularly monitored fields in January 1997, i.e., at the beginning of the second 
phase of the OGLE project (OGLE-II). The telescope was then equipped with 
the "first generation" camera with the SITe 2048 x 2048 CCD detector. The 
pixel size was 24 fj,m resulting in 0.417 arcsec/pixel scale. Observations of 
the LMC were performed in the "slow" reading mode of the CCD detector 
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with the gain 3.8 c~/ADU and readout noise of about 5.4 e~. Details of the 
instrumentation setup can be found in Udalski, Kubiak and Szymanski (1997). 

In June 2001 the third stage of the OGLE experiment (OGLE-III) began. 
The Warsaw telescope was equipped with a "second generation" CCD mosaic 
camera consisting of eight SlTe ST-002a CCD detectors with 2048 x 4096 pixels 
of 15 /Ltm size (Udalski 2003). This corresponds to 0.26 arcsec/pixel scale and 
the field of view of the whole mosaic 35' x 35'. The last observations presented 
in this paper were collected in May 2004. 

The vast majority of the observing points (430-860, depending on the field) 
were obtained through the /filter, while in the Vand i?-bands several dozen mea- 
surements were collected. The OGLE BV7 filters closely resemble the standard 
system. OGLE-II 7-band photometry was obtained using the Difference Image 
Analysis (DIA) method image subtraction algorithm developed by Alard and 
Lupton (1998) and Alard (2000), and implemented by Wozniak (2000). OGLE- 
III magnitudes come from the standard data pipeline (Udalski 2003). 1^ and 
5-band photometry was performed with the modified version of the DoPhot 
package (Schechter, Mateo and Saha 1993). OGLE-II and OGLE-III data were 
tied in the identical manner as in Soszynski et al. (2004). 

3. Selection of the Ellipsoidal Candidates 

We performed the period analysis for every star brighter than 7=18 mag. We 
used /-band light curves, for which the majority of observations had been col- 
lected and ran Fnpeaks program (Kolaczkowski 2003, private communication) 
to derive the most significant periodicities. 

For further analysis we selected stars with significant periodic light varia- 
tions, i.e., objects for which the ratio of the highest peak in the periodogram to 
the mean value was higher than 8. For this sample we constructed the logP-W/ 
diagram, where Wj is reddening free Wesenheit index, defined as: 

Wi = I-1.55{V-I) 

where /and V^are intensity mean magnitudes and 1.55 is the mean ratio of total- 
to-selective absorption {Aj / E{V — I)). It appears that logP-Wj diagrams can 
successfully replace period - near-infrared magnitude diagram, commonly used 
in the variable red giants studies. 

Our log P~Wi diagram shows a series of sequences, presented for the first 
time by Wood et al. (1999), and marked by letters A-E. In this paper we 
analyze the sequence E, spreading below the ridges of OSARGs, SRV, Miras and 
Long Secondary Periods (LSP), and extending down to our limiting luminosity 
(18 mag). 

We selected and visually inspected the light curves of stars forming the 
sequence E. It appears that most of these objects are characterized by nearly 
sinusoidal light curves with an amplitude between a few hundredths and a few 
tenths of magnitude. Fig. 1 presents exemplary light curves of these objects. It 
also appears that the majority of the light curves shows alternately shallower and 
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Fig. 1. Light curves of the exemplajry ellipsoidal red giants in the LMC. 
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Fig. 2. Light curves of the exemplary eclipsing red giants in the LMC. 

deeper minima, so the formal period is two times longer than the one obtained 
automatically in our period search. 

Among selected objects we also found 121 eclipsing variables, most often 
showing clear ellipsoidal effect. Several typical light curves showing eclipses are 
presented in Fig. 2. We compared the eclipsing binaries with the sample of 
sinusoidal variables, and found that when neglecting eclipses, the shapes of the 
light curves from both samples are very similar. 

Taking into account this similarity of the eclipsing and "sinusoidal" light 
curves, as well as the same position in the color-luminosity and period-luminosity 
diagrams, we interpret the "sinusoidal" objects as binary systems with orbital 
inclinations too small to cause the eclipses. Apparent variability of these objects 
is caused by the ellipsoidal distortion of red giants. 

To confirm whether a star is or is not an ellipsoidal variable, one needs to ob- 
tain a radial velocity curve. We believe that the variability of the vast majority 
of selected stars are caused by ellipsoidal modulation, but one should remem- 
ber that our sample may also contain a number of other types of variables. We 
expect such contamination especially among stars with the smallest amplitudes. 

It is worth mentioning that ellipsoidal light curves which, on first sight, 
do not show eclipses, can actually be eclipsing variables. First, there might 
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be grazing eclipses, which are not detectable by visual inspection. Second, 
even complete eclipses can be very shallow when the radius of the secondary 
component is significantly smaller than the radius of the red giant. Some of 
these objects can be X-ray eclipsing systems, which do not show eclipses in 
visual pass-bands, but they are observable in X-rays. 

4. Ellipsoidal and Eclipsing Variables 

We selected 1546 candidates for the ellipsoidal and 121 eclipsing variables in 
the LMC, with periods ranging from 15 to 600 days. Table 1 lists first 70 can- 
didates for the ellipsoidal variables from the field LMC_SC1. The following 
columns present star ID, star number (consistent with the LMC photometric 
maps, Udalski et al. 2000), equatorial coordinates, RA and DEC, for the epoch 
2000.0, periods in days (referred to the double-peak light curves), moment of the 
zero phase (corresponding to the deeper minimum), intensity mean IVB mag- 
nitudes, larger and smaller amplitudes of the /-band variability and remarks. 

The list of all variables is available in the electronic form from the OGLE 
Internet archive: 

'http:/ / ogle, astrouw. edu.pl/' 
ftp:/ /ftp, astrouw. edu.pl/ ogle/ogle2/var-stars/lmc/ ell^ 

or its US mirror 

http://bulge.princeton.edu/~ogle/' 
ftp://hulge.princeton.edu/ogle/ogle2/var_stars/lmc/ell/\ 

Individual 5 V7 measurements of all objects and finding charts are also avail- 
able from the OGLE Internet archive. The lists contain together 1753 entries 
but only 1667 objects, because 86 stars were detected twice - in the overlapping 
regions of adjacent fields. We decided not to remove twice-detected stars from 
the final list, because their measurements are independent in both fields and 
can be used for testing quality of the data and the completeness of the sample. 

In Fig. 3 we plotted the \ogP-Wi diagrams for the ellipsoidal binaries. 
The upper panel shows variables with the largest amplitudes (A/ > 0.08 mag), 
the middle panel presents variables with intermediate amplitudes 
(0.04 < A/ < 0.08 mag), while the lower panel - stars with the smallest am- 
plitudes (A/ < 0.04 mag). The dashed lines show the theoretical PL relations 
obtained assuming that the red giants fill up the Roche lobe. The mass of the 
red giant and mass of the companion were both set equal to IMq. In our cal- 
culations we used Paczynski's (1971) relation between the mean radius of the 
Roche lobe and the separation between components, which with the use of the 
Kepler law, leads to the following relation between the orbital period and mean 
density of the Roche lobe filling primary 



Porb = 0.12./5/(<z) [days] 
V P 
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Fig. 3. Period— VF/ diagrams for giants classified as ellipsodal variables. Upper panel shows 
variables with the largest amplitudes (A/ > 0.08 mag), middle panel presents variables with 
intermediate amplitudes (0.04 < A/ < 0.08 mag), while lower panel stars with the smallest 
amplitudes (A/ < 0.04 mag). 



where 

fiq) = V(l + g)(0.38-0.21ogq)3 

and q is the secondary to primary mass ratio. The /(q) dependence is weak: 
/ = 3 and 2.1 at g = 1 and 0.1, respectively. The values of Wi were calculated 
from Girardi et al. (2002) evolutionary tracks for the red giants phase at the 
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LMC composition (X = 0.74, Z = 0.008). Wc used Kurucz (1998) stellar atmo- 
sphere models to convert the evolutionary model parameters to the Wj values. 

It is clearly seen that dispersion of the sequence E depends on the amplitude 
of variability. The amplitude of ellipsoidal variation is a function of the third 
power of the Roche-lobe filling factor (Hall 1990), therefore the variables with 
the largest amplitudes should be close to the limiting lines. Periods of variables 
with smaller amplitudes of ellipsoidal modulation may significantly depart from 
the theoretical limits, because giants do not entirely fill the Roche lobe. 

The dispersion of data seen in the log P-W/ diagrams may be caused in part 
by the spread in the red giant masses. For instance, the models with M=1.5 Mq 
have the Wj values brighter by some 0.3 mag than the models with M = 1 Mq 
of the same Porb- Another contribution to the dispersion may arise from the 
the expected spread in the secondary mass, through the f{q) factor in Eq. (1). 

Different depths of minima, shown by the vast majority of our ellipsoidal 
candidates, can be explained by higer order effects of the tidal distortion. In 
addition, there is a number of variables in which different height of the light 
curve maxima is visible. This, so called O'Connel effect (Davidge and Milone 
1984), is usually explained as a result of spots in the stellar surface. In some 
cases maximum luminosity of these objects clearly changes with time, which 
is probably caused by the star-spot changes. Some unequal maxima may be 
a result of eccentric orbit of the system (see Section 6). 

5. Ellipsoidal Variability and Long Secondary Pe- 
riods 

Long Secondary Period (LSP) is one of the last unexplained types of stellar 
variability. At least 30% of the variable red giants show periods which are 
an order of magnitude longer than typical pulsation period of a semi-regular 
variable. Wood et al. (1999) noticed that the LSPs form a separate sequence in 
the period-luminosity space (sequence "D" ) . 

Many hypotheses have been proposed to explain the LSP. Radial and non- 
radial pulsation, periodic dust ejection, rotation of the spotted star and binarity 
have been suggested as the origin of the LSPs. Recently Wood, Olivier and 
Kawaler (2004) ruled out most of these propositions. They argued that the most 
likely explanation of the LSP is a low degree g+ mode of pulsation trapped in 
the outer layers above the convective envelope. 

We selected from OGLE database giants with dominant period correspond- 
ing to the LSP, and then we visually inspected their light curves. We removed 
from the list all uncertain and doubtful variables, and left about 1600 stars in 
which amplitudes of the LSP are significantly larger than amplitudes of other 
variations. Several typical light curves folded with the long periods are shown 
in Fig. 4. 

The period-W7 diagram for the dominant periods of these objects is shown 
in Fig. 5. Selected LSPs are marked with the empty circles, the magenta crosses 
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Table 1 

Candidates for ellipsoidal variables from the field LMC_SC1 (first 70 objects) 
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Remarks: U — uncertain, E — eccentric orbit, LMC_SC? — the same stars in the field LMC_SC? 
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Fig. 4. Light curves of the exemplary LSP red giants in the LMC. 



indicate our candidates for the ellipsoidal variables. It is worth noticing that 
apart from the P—L sequence D presented for the first time by Wood et al. 
(1999), Fig. 5 reveals additional sequence (marked by D') of stars with longer 
periods. 

It is not surprising that the ellipsoidal variability leads to a rather narrow 
band in the log P- 14^/ diagram. Along the red giant branch stellar radius, R 
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increases by more than two orders of magnitude. The luminosity increase is 
dominated by the radius increase and the period, which cannot be much longer 
than Porb given in Eq. (1), is determined primarily by the value of R^-^. 

The striking feature clearly seen in Fig. 5 is that the sequence of ellipsoidal 
variables partly overlaps with the sequence D and is a continuation of this 
sequence toward fainter stars. It might suggest binary origin of the LSP phe- 
nomenon in the red giants. 

Moreover, there is no doubt that some of the LSP cases may be explained by 
ellipsoidal variability. Fig. 6 presents exemplary light curves of stars where the 
characteristic OSARG variability is superimposed on the evident ellipsoidal light 
curves. We found about 300 such objects, mostly among brighter stars. The 
shorter, small amplitude variability corresponds to OSARG's period-luminosity 
sequences (Soszyhski et al. 2004). 

However, there are several serious differences between ellipsoidal variables 
and typical LSP giants. First, typical LSP light curve is not double humped, like 
the ellipsoidal light modulation (c/. Fig. 1 and Fig. 4). Second, the typical period 
- amplitude relation for the LSP variables is different than for the ellipsoidal 
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Fig. 6. Light curves of the eUipsoidal candidates exhibiting simultaneously OSARG-type vari- 
abiUty. Left column shows folded ellipsoidal light curves, and right column presents OSARG 
light curves after subtracting the dominant variability. 



giants. The amplitudes of LSP variability are positively correlated with the 

brightness of the star, while in ellipsoidal variables the amplitudes are connected 
only with the dispersion of the PL relation (Fig. 3). Third, Wood, Olivier and 
Kawaler (2004) noticed that the binary origin of the LSP is difficult to accept 
on statistical ground. LSP occurs in one third of the AGB stars, while the 
ellipsoidal variability relates to less than 1% of the analysed giants. 

The \ogP-Wi correlation for the LSP is very significant and cannot be ig- 
nored in any explanation of phenomenon. If the stars arc forced to rotate at 
nearly break up velocity, then we expect a -Prot(p) relation, the same as that 
given in Eq. (1), but with / = 1.2. Any asymmetry about the rotation axis could 
produce light variation with the Prot- However, the implied equatorial velocities 
of rotation seems unacceptably large. For instance at logP«2.5 we obtain the 
value of about 30 km/s. 

The explanation favored by Wood, Olivier and Kawaler (2004) (an excita- 
tion of a g-mode trapped above the convective envelope) may also explain the 
observed correlation. However, properties of such modes and chances of their 
driving remain to be studied. 
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Fig. 7. Light curves of the exemplary eUipsoidal red giants with eccentric orbits. 



6. Ellipsoidal Systems with Eccentric Orbits 

About 10% of candidates for ellipsoidal variables show evident deviations from 
symmetric sinusoidal light curves. For a sample of 165 objects the shapes of 
the light curves are completely different than typical sine-wave light curves of 
the ellipsoidal variables. Examples of these light curves are presented in Fig. 7. 
One should notice a large variety of the light curves. 

An explanation of such variables was proposed by B. Paczyhski (2004, pri- 
vate communication) who suggested that these objects may be ellipsoidal sys- 
tems with significant eccentricity of the orbits. A similar periodic light curve, 
although with the period of about 4 days, is shown by the well-studied star 
a Vir (Spica). Photometric and spectroscopic measurements showed that it is 
an ellipsoidal binary system with significant eccentricity of the orbit, equal to 
0.14. Variable distance of the binary components changes tidal interaction of 
stars during the orbital movement, which affects the shape of the ellipsoidal 
modulation. Moreover, the a Vir light curve is changing in the time scale of 
several dozen years (c/. the light curves obtained by Shobbrook et al. 1969 and 
Sterken, Jcrzykicwicz and Manfroid 1986), what is an effect of apsidal line ro- 
tation. Large variety of the light curves observed in the LMC can be explained 
by different longitude of the line of apsides. 
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Fig. 8. Model light curves of the ellipsoidal red giants with eccentric orbits. In the left 
column light curves with increasing eccentricity (e) are showed, the right column presents a 
series of light curves produced by binary systems with the eccentricity e = 0.2, and a different 
periastron length (p). 

To make sure that the variables are binary systems with eccentric orbits, we 
modeled such light curves using program Nightfall^ by R. Wichmann. A series 
of models are presented in Fig. 8. We modeled an equal- mass binary system with 
one of the components entirely filling the Roche lobe. The inclination angle is 
equal to 60°. In the left column the light curves with increasing eccentricity (e) 
are showed, while in the right column we present light curves produced by 
binaries with the same eccentricity {e — 0.2), but with a different periastron 
length (p). One can note clear similarity of the observed and calculated light 
curves. 

The models confirm that the larger eccentricity of the orbit, the more asym- 
metric light curve is. The variable showing one of the most asymmetric light 
curve - OGLE052013. 51-692253. 2 - corresponds to a model with e = 0.4. How- 
ever, precise modeling of the binary systems will be possible after obtaining the 
radial velocity curves. 

Verbunt and Phinney (1995) estimated that the circularization time scale of 
the semi-detached binaries containing a giant is of the order of several thousand 

\http://www.lsw.uni-heidelberg.de/users/rwichman/Nightfall.html\ 
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years. This estimation is almost independent of the mass of the companions. 
Thus, it is very unUkely that the close binaries maintained their initial eccen- 
tricity. It is more probable that the eccentricity of the orbits was induced by 
the rapid mass transfer or interaction of the third companion. 

Mechanisms of the orbit circularization induced by the tidal effects were 
studied by Zahn (1977, 1989) and Tassoul and Tassoul (1992). Such processes 
are especially effective for the red giants, due to the turbulent viscosity associ- 
ated to the convective envelopes. Studies of the spectroscopic binary systems 
containing red giants showed that there exist limiting period below which cir- 
cularization takes place. Mayor and Mcrmilliod (1984) analyzed the orbital 
parameters of 17 spectroscopic binaries with a red giant, and found that orbits 
are circularized at periods P< 150 days, however there are many exceptions to 
this rule. 

Obviously, the circularization time scale depends not only on the semi-major 
axis of the orbit (and consequently orbital period), but also on the radius of the 
giant star (related to the luminosity of the object). One can notice this effect 
in Fig. 5, where the period-luminosity relation of the candidates for eccentric 
ellipsoidal variables is marked by cyan dots. The magnitudes of objects are 
clearly correlated with the periods: the brighter stars, the longer orbital period 
of the system. Moreover, eccentric ellipsoidal variables have typically longer 
periods than non-eccentric ellipsoidal systems. 

7. Summary 

In this paper we showed that the ellipsoidal red giants constitute a numerous 

and homogeneous sample. This class of variables is difficult to detect because 
of their small amplitudes and long periods. Only long-term observing surveys, 
obtaining good quality photometry, can trace this type of variability. 

Ellipsoidal variables form the PL scqiience what is a projection of the radius - 
luminosity relationship for the red giants. We noticed that PL sequence of 
the ellipsoidal red giants is a continuation of the sequence of LSP variables, 
what may suggest a connection of the LSP and the binarity. Some ellipsoidal 
variables of the longer periods exhibit small-amplitude variability, corresponding 
to OSARGs. 

About 10% of the ellipsoidal variables show characteristic asymmetric light 
curves, probably caused by close binary systems with eccentric orbits. This 
sample may become an important test of tidal circularization theories. 
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